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Abstract 

Exploiting Virasoro constraints on the effective finite-volume partition function, we derive gen- 
eralized Leutwyler-Smilga spectral sum rules of the Dirac operator to high order. By introducing 
N v fermion species of equal masses, we next use the Virasoro constraints to compute two (low-mass 
and large-mass) expansions of the partially quenched chiral condensate through the replica method 
of letting N v — > 0. The low-mass expansion can only be pushed to a certain finite order due to de 
Wit-'t Hooft poles, but the large-mass expansion can be carried through to arbitrarily high order. 
Results agree exactly with earlier results obtained through both Random Matrix Theory and the 
supersymmctric method. 
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1 Introduction 



The exact spectral sum rules for the Dirac operator in finite-volume gauge theories derived by 
Leutwyler and Smilga ]|] and later by Smilga and Verbaarschot Q have sparked a huge activity 
in the field. It has turned out that there is an exact relation between universality classes in Ran- 
dom Matrix Theories and the three main ways in which chiral symmetry can break spontaneously, 
depending on the gauge group and the representation of the fermion fields ]3|, 0, H. The Random 
Matrix Theory formulation brings a completely different set of techniques into play, and much has 
been learned from this approach. One of the main efforts recently has nevertheless been to see how all 
the results obtained from Random Matrix Theory can be derived directly from the effective partition 
function. One of the first steps in this direction was a set of very compact relations that expressed 
microscopic spectral correlators (including the microscopic spectral density itself) directly in terms of 
effective partition functions with additional quark species ||. More recently, many of these results 
have been derived directly from the effective field theory [Q]. The idea has been to compute the (par- 
tially) quenched chiral condensate and higher chiral susceptibilities through an appropriately extended 
effective Lagrangian, and from this derive the microscopic spectral correlators. As the technique of 
ref. |?J makes use of additional quark species, of which half are bosonic and the others are fermionic, 
this extends the flavor symmetry to a super Lie group at intermediate steps. For this reason it is 
commonly known as the supersymmetric method (although it, as applied, has nothing to do with 
space-time super symmetry). 

The purpose of this paper is to explore another technique that can be used to derive the partially 
quenched chiral condensate. This is based on the so-called replica method^ in which one extends the 
theory with N v additional fermionic species, and takes the limit N v —* at the end of the calculation. 
The replica method is known in other contexts to be often problematical [|| (see however also the 
recent discussion on this issue ||), but we shall find no fundamental difficulties with this technique in 
the present context. 

Let us first outline how the replica method can be used to compute partially quenched averages in 
QCD. Consider the QCD partition function with Nf physical quark fields and N v additional quark 
fields of degenerate masses m v . We restrict ourselves to gauge field sectors of fixed topological charge 
v, which we for simplicity from now on take to be non-negative: 

(Nf \ N f 
[] m v f m£ vV f[dA] u det'(# - m/" J] det W ~ m f) e" s ™^ , (1.1) 
/=i / /=i 

where the determinants are taken over non-zero modes only. This can be considered as an average 
over gauge fields (and physical fermions, here already integrated out) of N v identical replicas of the 
fermionic partition function 



Z v = J dipdip exp J d^xi^iil/) — m^ip 



(1.2) 



in the sense that 



/ N f \ N f 

z (N f+ N v ) = JJ m , J [dA]v "Q detW _ mf) [Zv] N v e -S YM (L3) 



Or replica "trick" , a name that suggests magic and trickery. We prefer the more neutral terminology. 
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In condensed matter physics this (unnormalized) average is conventionally denoted by a bar: 



(1.4) 



The subscript v will throughout denote "valence", and the additional quarks are thus valence quarks, 
a terminology borrowed from lattice gauge theory. Obviously, if we let N v = we simply recover the 
original QCD partition function. For Nf + N v not too large, the theory is presumed to undergo, in the 
chiral limit of massless quarks, spontaneous chiral symmetry breaking according to the usual pattern 
S\J(N f + N v )xS\J(N f + N v ) -> S\J(N f + N v ). Taking the limit N v -> must therefore incorporate 
some kind of analytic continuation in re of the Lie group SU(n). 



For any fixed N v we can compute observables in the extended theory by viewing the partition function 
as the generating function of the re-point functions of ipip. We shall here focus on just the chiral 
condensate itself, and only on the chiral condensate of the additional fermionic copies. As we will take 
N v — > 0, this becomes the partially quenched condensate in the theory with Nf physical fermions. We 
define this partially quenched (mass-dependent) chiral condensate by 



MhlMl s lim J-^-ln^^ 



(1.5) 



where S is the physical infinite- volume chiral condensate in the theory with Nf fermions, fij = rrijVTi 
and similarly fi v = m v VTi. Treating N v as a parameter that is not restricted to be integer (which is 



permitted by the representation of the partition function in eq. (1.1)), we can make a Taylor expansion 
in N v : 



~(N f +N v ) 
and thus write, equivalentlyj^j 



z, 



(N f ) 



+ N v —Z; 

oN v 



d (N f +N v ) 



N v =0 



>(*» 



- 1 _d d_ (N f +N v ) 

dN v d\i v 



N v =0 



(1.6) 



(1.7) 



The partially quenched chiral condensate is a particularly convenient quantity for lattice gauge theory 
simulations, and comparisons with theory have already been made both in the topologically trivial 
sector [10] and, very recently, in sectors of non- vanishing topological charge v [JO]]. We note here 
that higher-order partially quenched chiral susceptibilities, which are also readily studied by Monte 
Carlo simulations, can be derived in a completely analogous manner. For a chiral Appoint function one 
simply needs to introduce N v \, . . . ,N v k sets of different additional species (each of iV„j'th degenerate 
masses ftm), and take the combined limit of all N v i — * in the end, after having performed the required 
differentiations. 



In this paper we will show how the replica method can be applied to the effective finite- volume partition 
function, and in this way we shall derive analytical series expansions for the partially quenched chiral 
condensate. As a by-product of the analysis we will get, for free, generalized Leutwyler-Smilga sum 
rules to very high order. The plan of the paper is as follows. In the next section we review the 
technique for deriving high-order expansions of the effective partition function on the basis of an 
iterative sequence of partition function constraints. From this expansion we immediately derive a long 
list of generalized spectral sum rules of the Dirac operator. In section 3 we turn to the replica method, 
where we first exploit the same partition function constraints to derive a small-mass expansion for the 

2 One can explicitly check that the two derivatives commute in our case. 
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partially quenched chiral condensate. We also show how a different set of partition function constraints 
can be used to derive a very high order expansion in large masses. In section 4 we point out that the 
usual Leutwyler-Smilga spectral sum rules, which by definition are taken with respect to the massless 
theory, have a natural generalization where the spectral sums are taken with respect to the theory with 
massive fermions. We derive a series of such massive spectral sum rules, perturbatively expanded in 
the physical masses. In section 5 we discuss some generalizations to different patterns of spontaneous 
chiral symmetry breaking, and point out why the present series expansions cannot be used to derive 
the microscopic spectral density of the Dirac operator itself. In section 6 we briefly list the main 
conclusions. 



2 Spectral Sum Rules from Virasoro Constraints 



Spectral sum rules of the Dirac operator are derived by comparing the full QCD partition function 
with the effective low-energy partition function. In the large- volume scaling region V <C this 
effective partition function is to leading order JIJ 



(N f ) 



ueU(N f ) 



dU (det U) y exp 



\tt(MU^ + UM ] ) 



(2.1) 



where M is the quark mass matrix rescaled by the space-time volume V times the value £ of the 
chiral condensate in the chiral limit (and at 6 = 0). We will in fact always take this to be diagonal 
of entries /ij = rriiYiV. In this section we shall restrict ourselves to Nf (physical) fermions with no 
additional quark species, and N v will thus here be taken to be zero from the very beginning. The 
integration in eq. ( |2.1| ) is over the coset of chiral symmetry breaking, here extended from SU(iVj) to 
\J(Nf) due to the projection on a sector of fixed topological charge v. The integral is known in closed 
form p|], 



det 
A({/. 2 }) 



where the Nf x Nf matrix A is given by 



(2.2) 



(2.3) 



I n {x) is a modified Bessel function, and the denominator is given by the Vandermonde determinant 
of rescaled masses: 



N f 



A({/u 2 }) 



n^ 2 

i>3 



det 

i,3 



?y-i 



(2.4) 



Although the effective partition function is known explicitly, eq. (^^) is not in a form suitable for the 
derivation of spectral sum rules of the Dirac operator. To get it into such a form, it is convenient to 
start with the case v = 0. In that case the left and right invariance of the Haar measure under unitary 
transformations shows that the partition function depends only on the combination M^M. A crucial 



observation of ref. [13, 14] is that the partition function in fact satisfies an infinite set of constraint 
equations that can be used to determine it uniquely. The precise form of these constraints depend on 
the chosen variables, and one can use two different sets: 

t+ = ±-Tr(M^M) k , k = l,2,... , (2.5) 
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which is a suitable set for small-mass expansions, and 

2 2fc+l 

£ = -^ TI Tr((A<t A< )-(2fe+i)/2) j = ,l,... , (2.6) 

which is a suitable set for large-mass expansions (the coefficients in front have been inserted for later 
convenience). Spectral sum rules of the Dirac operator are derived by means of a small-mass expansion, 
and we shall therefore begin with the description in terms of the variables^ ft . 



Defining, for n > 1, 

o oo o n— 1 02 

OTn fc=0 OT n+fc fe=l OT fc OT n-fc 

the partition function ( p. 10 ) is found to satisfy [14| 

C+4" f) = *nA Nf) ■ (2-8) 
These constraints are consistent in the sense that they satisfy the classical Virasoro algebra 

[£+,£+] = (n-m)Ct +m , (2.9) 

so that they do not generate new constraints beyond the infinite tower of Ct 's. The constraints are 
also complete: They determine the partition function uniquely, given the boundary condition that 
Z^ Nf) = 1 for all 4 = 0. 

(N f ) 



The small-mass power series for Z can conveniently be chosen [14] 

4^ = i+E E ^(W) r rl'"'+t^ ' (2 - 10) 

where use has been made of the boundary condition to determine the zeroth order coefficient. 

The coefficients C^A{k\) are determined uniquely by the Virasoro constraints They factorize 

into a polynomial part, C^A{k\), and a singular part []T4| : 

K({fc})-1 

({*;}) = ({*:}) LI (^/-^T 1 , (2-11) 

where K({k}) = J^i.ki- The singular parts will play an important role in what follows. Originally 
noted by de Wit and 't Hooft in the context of strong-coupling expansions in lattice gauge theory |l5| , 
these poles are in fact innocuous as far as the expansion of the partition function itself is concerned 
pf. What happens is that for (integer and non- vanishing!) values of N the poles only occur in 
terms that are linearly dependent, and therefore should be combined. After combining them, all poles 
cancel. This phenomenon becomes obvious when we consider the explicit expansion below. 



The coefficients CN f ({k}) have already been computed recursively in ref. |14j to high order, enough 
to determine the partition function to 8 th order in the masses. The restriction to the v = sector can 
easily be lifted due to flavor-topology duality, which states that 

zi Nf \M,M*) = [det(M)] u Z { Nf+u) (MM^) (2.12) 



3 Since M f M is a hermitian matrix of size Nf x Nf, only Nf of these variables are independent; this is of no 
consequence for the subsequent analysis. 



5 



where Zq N * +u \mM^) is the v = partition function extended with v massless quarks [18]. Intro- 
ducing N = Nj + N v + v (here with N v = 0) we have 

1 ._,..+ . •> 1 



2f'> = det-(M) (l + W ^M + __ _ _ {T r(M*M)) 2 - 3m{N2 _ ^ Tr^M) 2 

N 2 — 2 5 

-(Tr(.M t .M)) 3 ^ ^ rTr((M*M) A ) 



384N(N 2 - l)(N 2 - 4) v v " 1024iV(iV 2 - \){N 2 - 4)(iV 2 - 9) 
IN 2 — 3 

-T¥(7W t 7W)TV((7W t 7W) 3 ) 



768iV 2 (iV 2 - l)(iV 2 - 4)(iV~ 2 - 9) 



+ N2 + 6 (Tr(M^M) 2 ) 2 

^2M8N 2 (N 2 - 1)(N 2 - 4)(iV 2 - 9) { { ' ' 



1024iV(iV2 - l)(iV 2 - 9) 

+ *imw^w^^ mMiM)),+ --j <2 ' 13) 

We have checked that this expansion is correct to the order given by expanding, for fixed Nf, the closed 
expression fl2.2|) up to that order. (Note that it is not at all simple to rearrange the closed expression 
( |2.2j ) in terms of an expansion in just powers of traces of the mass matrix; the rearrangements are 
different for each value of Nf). In this expansion we directly verify the cancellation of de Wit-'t Hooft 
poles for any of the finite, integer, values involved. 

With this high-order expansion of the effective partition function, we can now easily derive a long 
series of spectral sum rules, following the method of Leutwyler and Smilga Q. Because the method is 
described in detail in ref. jjj, we shall not give details of how these sum rules are extracted, but only 
quote the results. We rescale as usual the eigenvalues Aj of the Dirac operator according to Q = AjEV. 
Expanding the original partition function ( |1.1| ) around the massless theory, and comparing with the 
above expansion up to 8 th order, we then find: 

ih <2 - 14) 

(2.15) 




16 iV(iV 2 - 1) 
1 1 



32 N(N 2 - 1)(N 2 - 4) 

5 1 
256 N(N 2 - 1)(N 2 - 4)(iV 2 - 9) ' 



(2.16) 
(2.17) 



E c2 

Sri 



vCn>0 



E 

V Cn>0 



c 2 



1 



1 



16 N 2 - 1 



1 N 2 — 2 

MN(N 2 - l)(N 2 - 4) 



(2.18) 
(2.19) 
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N 2 + 6 



256 N 2 (N 2 _ i)(AT2 _ 4)(at2 _ 9) 



(2.20) 



E 



1 



c 2 C A 

1 



E 2 , (i 

\Cm,Cn>0 ^rnSn, 



E 

\Ck,tm,Cn>0 



E 

^Cfc>Ci,Cm>Cn>0 



1 



< A2 r 2 r 4 

Sfc SmSn 



/•2A2A2 Z-2 



64 (N 2 - 1)(N 2 -4) 

1 2iV 2 - 3 

256 iV 2 (iV 2 - l)(iV 2 - 4)(iV 2 - 9) 

1 1 

256 N(N 2 _ i)(AT2 _g) 

1 N 4 - 8N 2 + 6 

256 N 2 (N 2 - l)(iV 2 - 4)(iV 2 - 9) 



(2.21) 
(2.22) 
(2.23) 
(2.24) 



Of these, the relations ( 2.14 ), ( 2.15| ) and ( p. 18 ) have been derived previously Q (see also ref. [I7j] ). We 
note that the DeWit-'t Hooft poles make these sum rules divergent for an increasingly larger number 
of integer iV- values. This is due to infrared divergences near £ ~ 0. Just as these poles cancel among 
different terms in the small-mass expansion, one can form analogous combinations of sum rules that 
are infrared finite. For example, the de Wit-'t Hooft poles at N = 1 associated with the terms of 
order /i 4 in eq. ( |2.13 ) cancel by combining the third and fourth terms in that expansion. This directly 
translates into a cancellation between two eigenvalue sum rules from our list: 




1 



16JV(JV + 1) 



(2.25) 



a cancellation that was already noticed by Leutwyler and Smilga Similarly, while sum rules (2.16), 
( P19D and ( ggj] ) are individually singular at both N = 1 and N = 2, the combination 



\Cn>0W \\Cn>0^ 



E 



1 



1 



\Cm,Cn>0 



c 2 C A 



64N(N + 1)(N + 2) 



(2.26) 



is finite (this is 6 times the 0(1/ A 6 ) contribution to the partition function). A similar phenomenon 
persists to all orders. 



The massless spectral sum rules (|2~ll)-(|2~24|) can of course be compared with the chiral Random 
Matrix Theory result. This is done by averaging the inverse moments according to the universal 
massless microscopic spectral density Q 



^"'(0 = I 1 r 



JN f +v((Y 



e.g. 



E 

\(n>0 



2A- 



■^v / +i/-i(C)^jv / +i/-i-i(C) 



°° At gg(0 
-00 S 



(2.27) 



(2.28) 



We have explicitly checked a number of the above sum rules, always finding completely agreement 
with the Random Matrix Theory results. 
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3 The partially quenched chiral condensate 



Let us now extend the effective theory with N v valence quark fields. We are interested in obtaining the 
partially quenched chiral condensate, not higher order susceptibilities, and therefore take the valence 
quark masses as degenerate and of size fi v = m v T,V. As mentioned above, the extended partition 
function is also the generating functional for the n-point functions of ipip (quenched or unquenched); 
the partially quenched chiral condensate is obtained through the replica method. Since the partition 
function for v = only depends on Ai^Ai and the valence quarks have degenerate masses it is natural 
to introduce a mass matrix of the physical quarks M so that 

Tt(M^M) = Tr(M 2 ) + N v n 2 v . (3.29) 

In short, for the definition of the partially quenched chiral condensate we can take 

^ lim J-^lnzfW . (3.30) 

The subscript on T, refers to the topological sector, as usual. As we have seen above, it is not 
particularly helpful to know the exact analytical expression for the (extended) effective partition 

function Zu Nf+Nv ^ if one wishes to derive spectral sum rules for the Dirac operator. There a small- 
mass expansion is required by definition. As for finding the partially quenched chiral condensate using 
the expression ( 3.30) ), no expansion in either large or small masses is required in principle. However, 



the closed analytical formula ( |2.2| ) is unfortunately not directly suitable from our point of view. This 
is because the ^-dependence enters in a quite non-trivial way through the size of the matrix whose 
determinant needs to be taken. Because there is no simple extension of the analytical expression (2.2) 



outside integer values of N v , we restrict ourselves here to series expansions. We shall in fact be able 
to carry such expansions through in both the large-mass and small-mass regions. 

3.1 The small-mass expansion 

We begin with the small-mass expansion, because here we already have the needed expansion at hand 



see eq. ( 2.13 )). Our first observation is that N v enters in a manner that allows for an analytic 



continuation once we make use of eq. fl3,29| ). We can thus proceed with the replica method. 



Inserting the expansion fl2.13|) into ( 3.30| ) we obtain the partially quenched chiral condensate to 7 



th 



order in the masses (truncating the expansion at this order is a matter of choice; we do it because 
it corresponds to consistently expanding both in the valence quark mass [i v and the physical fermion 
masses fii to the same order): 



^ — + 



1 TrM 2 (TrM 2 ) 2 



2N 8N 2 (N 2 - 1) 8N 3 (N 2 - 1)(N 2 - 4) 

TrM 4 15TYM 6 
+ 



16N 2 (N 2 - 1)(N 2 -4) 384N 2 (N 2 -l)(N 2 -4)(N 2 -9) 

3(2N 2 - 3)(TrM 2 ) 3 3(3 - 2 N 2 ) TrM 2 TrM 4 

+ 32N 4 {N 2 - 1) 2 (N 2 - 4){N 2 - 9) + 32N*(N 2 - 1) 2 (N 2 - A)(N 2 - 9) + " ' 



-/4 



1 TrM 2 
+ 



8N(N 2 -1) 8N 2 (N 2 - l)(N 2 - 4) 
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+ 



(6iV 2 - 9)(TrM 



2^2 



-Hi 



32^3(^2 _ i)2(jv-a - 4)(iV 2 - 9) 

3(3iV 2 - 7)TrM 4 
12m 2 {N 2 - 1) 2 (N 2 - 4){N 2 - 9) 

1 



+ ... 



+ 



15TrM 2 



16iV(iV 2 - 1)(N 2 -4) 128iV 2 (iV 2 - l)(iV 2 - 4)(7V 2 - 9) 
5 



128N(N 2 - 1)(N 2 -4)(N 2 



+ 



(3.31) 



where N = Nf + v (since the N v — > limit has already been taken). Note that the "topological" term 
v / H v comes out trivially from the (det .M) v -factor in front, as it should. 



Is this expansion correct? There are by now two independent ways to check it: Comparing with 
Random Matrix Theory, and comparing with the supersymmetric method of partially quenched La- 
grangians. Before proceeding with such comparisons, let us however first note an obvious feature of 



the above expansion. In contrast to the small-mass expansion of the partition function itself ( 2.13 ), 
there is clearly not going to be a cancellation of the de Wit-'t Hooft poles in the expansion of the 
partially quenched chiral condensate. This may appear surprising, since the chiral condensate is just 
determined from the partition function through, for example, a formula like ( |1.5| ). If the de Wit-'t 
Hooft poles cancel in the partition function, how can they not cancel after taking the //^-derivative 

1.5)? The answer is simple. To be able to take the quenched limit 

This 



as in eq. 



and the limit N v 

of N v — > we must not cross any singularities that can obstruct the analytical continuation. 



is, however, precisely what happens for those terms in the expansion flOll) tha t, for given N, are 
singular. To give an example, let us consider the /i^-term in the expansion ( 3.31 ). For simplicity, let 
us also consider the case where all physical masses vanish. Then there is only one term: 



8N(N 2 - 1) 



(3.32) 



This term is singular for N = 1 (and, of course, also in the more trivial case of N = 0), and therefore 
cannot possibly represent the partially quenched chiral condensate to order ^ in the N = 1 theory. Let 
us therefore trace what happens if we set Nf + v = 1 from the outset. The only term in the expansion 
( f2. 13 ) that gives an 0(fj%) contribution to T, u (fi v ) in this massless case is (N = Nf + N v + v) 



32N(N 2 - 1 



-Tv(M^Mf 



32N(N 2 - 1) 



32(N v + l)(N v + 2) 



(3.33) 



which precisely in this Nf + v = 1 case is not linear in N v . The factor 1/N V in the expression for 
the quenched condensate ( 3.3C| ) will in this case remain uncancelled, and an analytic continuation to 
N v = is prohibited by the singularity. The term ( |3.32| ) in the expansion ( |3 .31 ) is therefore incorrect 
precisely in the case Nf + v = 1 (but it is valid for all Nf + v > 1). This holds also for all higher 
orders in the expansion. The small-mass expansion for Ti u is thus given to a finite number of terms. 
For example, in the massless Nf + v = 3 theory we obtain 



s 



fi v 6 



192 1920 "' 



Nf + v 



(3.34) 
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where the higher-order terms in fa cannot be purely power-like. The same phenomenon occurs for 
larger integer values of N, where we encounter the singularity at higher orders. We thus conclude that 
the small-mass expansion ( 3.3l[) is valid up to the first term where a pole appears. The obstruction 



to an all-order expansion of H v {fa, {fa}) is from this point of view uncancelled de Wit-'t Hooft poles. 
The expansion can of course be pushed to very high orders by considering Nf + v sufficiently large. We 
stress again that the expansion is perfectly meaningful up to the order at which the first uncancelled 
de Wit-'t Hooft pole appears, and that it is completely understood why the expansion cannot be 
pushed beyond this order. It is not a peculiar and unwanted artifact of the replica method, but a real 
feature of the small-mass expansion of the partially quenched chiral condensate. 



With this caveat in mind, we can now check the expansion ( 3.31| ) against both the supersymmetric 
method and results from Random Matrix Theory. For instance, both methods yield a fully quenched 
chiral condensate of the form ^] 

% v = fa %{fa)K v {fa) + I v+ i{fa)K v -x(fa)] H . (3.35) 

hi fa; 

where I n (x) and K n (x) are modified Bessel functions. Looking at the small-//,, expansions of K n (fa), 
one notices that they are only purely power-like up to order n — 2 (after which a logarithmic term must 
be included). The small- mass expansion of the closed expression ( |3.35D therefore cannot possibly match 
to all orders the purely power-like expansion that would seem to follow from the small-mass expansion 
of the partition function in eq. ( 2.13| ). And we find that the agreement of the two expansions is exact 



precisely up to the order at which our small-expansion (3.31) ceases to be valid. The appearance, at 
higher orders, of logarithmic terms in the expansion explicitly confirms the statement above: there 
is simply no power-law small-mass expansion of the chiral condensate beyond that given order. It is 
quite remarkable that the normally cancelling de Wit-'t Hooft poles precisely obstruct the small-mass 
expansion at just the right order. 

We have also checked that the double-expansions in both fa and physical fermion masses fa are correct 
precisely up to the order at which eq. (3.31) is valid. We have done this by comparing the small-mass 



expansions of the partially quenched chiral condensate in the Nf = 1 theory with the expansions of 
the closed expression 



^v{fa,fa) = [I u+l ^ v )K u+ i(fa) + Iv+2{fa)K p (fa)} + — 

h fa 

v) fajlv 

(fa)I u+ i(fa) - faT v (fa)I v+1 (fa 

i v (n) i4 - n 



which has been derived by the supersymmetric method (7j (and which can also easily been seen to 
agree with the Random Matrix Theory result). Again we find perfect agreement with all terms up to 
the order at which our expansion ( 3 . 3 1| ) ceases to be valid, i.e. up to the first uncancelled de Wit-'t 
Hooft pole. 

3.2 The large- mass expansion 

Having succeeded in deriving the low-mass expansion for the partially quenched chiral condensate, we 
next turn to the opposite expansion, i.e. for large masses. The advantage of this expansion is that it 
can be pushed to arbitrarily high order, for any value of Nf. 

Recall that the effective partition function ( |2.lD only depends on the combination M^M for v = 0. 

(JVf) 

This combination is hermitian and there are thus only Nf degrees of freedom in Z , represented 
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for instance by the eigenvalues /xf of M) M.. In terms of /i, the expansion variables of eq. ( |2.6| ) read 



22fc+l 



TV 



= Y^P 2k+L> . (3.37) 

The idea is now to find suitable partition function constraints that will enable us to solve for the 
partition function in an expansion in the variables tZ- This is not completely straightforward, but as 
shown first by Gross and Newman jy]] , it is possible to recover a complete set of Virasoro constraints 
in the variables tZ if one first extracts a simple prefactor from the partition function. The factorization 



is as follows JL3|] (see also ref. [14] for the generalization to other types of matrix integrals): 

Nf 

1 2' 



4 Nf) = (ml^ + Ur 1/2 )e^Y Q ^\{t-}) . (3.38) 



o,6 



The exponential prefactor indicates that an expansion of the partition function entirely in terms of the 
t^T 's simply is not possible. This is also obvious if we consider the original group integral in a saddle- 
point approximation around large A4: the leading behaviour is indeed exponentially large. Moreover, 
one suspects from the form ( p. 38 ) which is essentially a saddle-point expansion that the resulting series 



expansion we will derive below is at best asymptotic. This is in contrast to the previously considered 
small-mass expansion, which we believe is convergent. 

The remaining factor Yq ({tZ}) in eq. ( |3.38| ) turns out to be annihilated by an infinite set of 
Virasoro operators. Explicitly, by defining 



k =o - 9t k 16 dt o 



Y( k + L) t -^_ + Ly ° + n >i (3.40) 

^ 2 k r)i~ 4 ^ r)t~ r)t~ f)t~ 

k=0 OZ k+n * k=l ° l k-l m n-k ° ln 



1 n d 1 d , . 

- Y.(k+- )t k ^ + - + 7 ^, (3.39) 

the function Y^ Nf \{tZ}) of (gjD is found to satisfy @ 



C~Y^ Nf) =0 n > . (3.41) 

One readily verifies that these constraints indeed also fulfill the commutation relations of classical 
Virasoro generators (2.9), and thus form a consistent set of constraints. These constraints are also 



complete in that they determine the partition function uniquely, given the boundary condition Y = 
1 for tZ = 0, i.e. in the limit of infinite masses. While these constraint equations were already 
established in ref. [13], they have not previously been used to derive a systematic expansion for the 
partition function. We do this by copying the procedure of the small-mass expansion. That is, we 
expand the unknown function YQ Nf \{tZ}) as follows: 

oo 

Y iNf \{tZ}) = 1 + E C ^+ E <W^ fc "+... (3.42) 

fc=0 0<fci<fc 2 



We next solve for the coefficients c^. 1 ... i- n iteratively. In fact, this is in many ways simpler than in 
the small-mass expansion in that we in this case easily can derive closed analytical expressions to all 
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orders. For example, we find the following simple formula for the first strings of coefficients: 

(fr -1 8m(2k + l)+8j + l \ 

Co ,o_o, fc , m ., fc = yy[ -"t 16 J c ^ ■ ( 3 - 43 ) 

Using this formula, and others similar, it is straightforward to derive very high order expansions for 

^0 (te })' anc ^ therefore also for the partition function itself. For example, to 7 th order in the 
masses we find a rather formidable- looking expression (that is easily pushed to much higher orders): 

N f 



f 1 1 

4 N,) = (ntU+iU) 



, , , , " 1/2 )exp 



1 + iMM 1 MY 1 ' 2 + — (Tr(M*M)~ 1/2 ) 2 
8 128 

+ I | i (Tr(^t^)-i/2 ) 3 + ATr^t^)-^ + ^{^M)-^f 
+ T ^ i (Tr(Mt^)-i/2 )(Tr( ^t^ ) -3/ 2) + J|^(Tr(^t^)-i/2 )5 

+^(Tr(^t^)-i/2)2 (T¥( ^t^ ) -3/ 2) + -^Tr(M^M)-^ + 4^(Tr(^^)- 1/2 ) 6 
+^(^(^t^)-V2 ) 3 (Tr( _ M t A1) -3/ 2) + ^(Tr^t^)-3/ 2)2 
+^(Tr^t^)-i/2 )( T r( ^t^)-5/ 2) + J|^|(TY(^t A1) -i/2 ) 7 

+ 4iIii (Tr( - Mt - Mrl/2)4(Tr( - Mt - Mr3/2) + ^(^(•^ t ^) -1/2 )(^(^ t ^)" 3/2 ) 2 

+|^(Tr(^t^)-V2 ) 2 (Tr( ^t^ ) -5/ 2) + ^Tr(M^M)-^ + . 



(3.44) 



In writing fl2.6|) and ( 3.37 ) we have implicitly assumed that ftA^Ai is invertible. This means that 



there must be no zero-eigenvalues in the mass matrix, and we thus by default cannot consider any 
number of massless fermions in this large-mass expansion. In particular, the concept of flavor-topology 
duality, by which a topological charge v is implemented by adding v massless fermions to the v = 
partition function, cannot be applied in the framework of a large-mass expansion. Because the Virasoro 
constraints ( ET39D and ( gig ) apply only to the case v = we know of no analogous way to implement 
this large- mass expansion outside of the v = sector. 

Introducing N v mass-degenerate valence quarks into the large-mass expansion of the effective partition 
function, and noting that in this case 

Tr(M^M)- 1/2 = Tr(Af)- 1 + iV> _1 , (3.45) 

we again observe that N v can be regarded as a continuous parameter. The replica approach can 
therefore again be applied. 

Using the definition ( |3.30| ) and a good handful of simple algebra we find the partially quenched chiral 
condensate for large masses: 

goQ^M) = 1 1 
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P-v 



1 1 ,1, If 1 

8 + 8^M> + 8 ( M 



9 m / 1 V 

+ 128^M3) 



1 / , 1 A 3 1 / , 1 A 4 9 , 1 N ,1, If , 1 V 



45 m , 1 . / , 1 



+ 1024 Tr( M^ ) + --- 



1 

"/4 



/'(' 



9 



27 



1 



27 



64 



45 



128 + 128 Tr( M ) + ^( Tr( ^)) + 5ZR(lf)) +^ Tr W + 



6 -J 



189. 



1 



512 V M 3 



225 1125 m , 1 . 

1 Tr( — ) + . . . 

1024 1024 y M J 



(3.46) 



This expansion agrees with the asymptotic expansion of the analytical expressions for {^}) found 

in @ for N f = and N f = 1. 

4 Massive spectral sum rules 



The spectral sum rules presented in section || concern inverse moments of the Dirac eigenvalues aver- 
aged with respect to the massless theory. These massless spectral sum rules are, however, only special 
cases of massive spectral sum rules obtained by averaging with respect to the massive theory. In fact, 
massive spectral sum rules appear quite naturally in lattice simulations. Such massive spectral sum 
rules are conventionally taken H |l^] to be of the form 



1 



(C 2 + ^ 



for given integer re/s. Here masses and Dirac operator eigenvalues are treated on equal footing in the 
denominators. However, we could equally well define massive spectral sum rules by 




ps(C; {Pi}) 



C 



2n 



(4.47) 



which are just the usual inverse moments, but evaluated in the massive theory. Here and below the 
brackets denote the average with respect to the massive and partially quenched theory in a topological 
gauge field sector of charge v. In order to derive such massive spectral sum rules we return to the 



small-mass expansion of the partially quenched condensate ( 3.31 ). 



Inserting the original partition function ( |1.1| ) into the definition ( |3.30| ) of the partially quenched chiral 
condensate one has 



^u(Pv,{p}) 




(4.48) 



We can now find a useful general relation between the partially quenched chiral condensate, and 
(massive) spectral sum rules. The connection is simple: the sum rules of inverse moments £ _2n of the 
eigenvalues will be convergent up to a given value of n, which we denote by k. Let us therefore expand 
the denominator of eq. ( |4.48j) in partial fractions up to this maximal value k: 



1 



c 2 c 4 

Sn Sn 



+ 



(-1) 



D k pf 



Cn^iCn + Pv 



(4.49) 
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which means that 



s Vis, « / * Vf?o ci / * ' " Vfs, era + / * • (450) 

We can thus simply read off a whole string of massive spectral sum rules from the coefficients in the 



//,,-expansion of the partially quenched chiral condensate (3.31): 



' 1 \ 1 TrM 2 (TrM 2 ) 2 

^ nCl/ ~ 4iV + 16iV 2 (iV 2 - 1) + 16iV 3 (iV 2 - l)(iV 2 - 4) 

Tr((M 2 ) 2 ) 15Tr((M 2 ) 3 ) 



+ 



32iV 2 (iV 2 -l)(iV 2 -4) 768iV 2 (iV 2 - l)(iV 2 -4)(iV 2 
3(2iV 2 - 3)(TrM 2 ) 3 



64N 4 (N 2 - l) 2 (iV 2 - 4)(iV 2 - 9) 
3(3-27V 2 )TrM 2 Tr((M 2 ) 2 ) 
+ 64iV 3 (iV 2 - 1) 2 (7V 2 - 4)(iV 2 - 9) + " ' ( ' 51) 

' ^ 1 \ 1 TrM 2 

nCl/ ~ 16iV(iV 2 - 1) + 16iV 2 (iV 2 - 1)(N 2 - 4) 

(6iV 2 - 9) (TrM 2 ] 2 

\(N 2 - 9) 

+ . . . (4.52) 



64A^ 3 (iV 2 - l) 2 (iV 2 - 4)(iV 2 - 9) 
3(3iV 2 - 7)Tr((M 2 ) 2 ) 



256iV 2 (iV 2 - 1) 2 (7V 2 - 4)(iV 2 - 9) 
/ „ 1 \ 1 15TrM 2 

V nCl/ ~ 32iV(iV2 - 1)(AT2 _ 4) + 256iV 2 (iV 2 - l)(iV 2 _ 4)(iV 2 - 9) + " ' ( ' 53) 

(^ Q Cl) = 256iV(iV 2 -l)(iV 2 -4)(iV2_ 9) + --- (454) 

The resulting sum rules are here only given as a perturbative expansion in the physical fermion masses, 
in contrast to the massless spectral sum rules which are exact. One notices that trivially the "diagonal" 



massless spectral sum rules ( 2.14 )-( 2.17 ) are reproduced by taking all ^ = 0. As discussed already in 
section 2, the above expansions are valid for N larger than the biggest integer poles. This is completely 
understandable from the present point of view, as it corresponds to performing the expansion ( |4.49 ) 



in partial fractions only up to the point where the spectral sums of all terms still converge. The 
"non-diagonal" massive sum rules can also be calculated within this framework. One simply needs 
to break the mass-degeneracy of the valence quarks and consider general partially quenched n-point 
correlators of ipip, as discussed in the introduction. 



These new massive spectral sum rules can easily be compared with the predictions from Random 
T 

of TrM 2 , and then insert this expansion in 



Matrix Theory. What we need is to expand the microscopic spectral density (C; {/^i}) m powers 



^■w-^r^fni. (4.5 5 ) 



S Jo C + K 

For example, for Nf = 1 the microscopic spectral density for a massive fermion in a gauge field sector 



14 



of arbitrary topological index v is given by (T{|| (see also the second reference of |g] and ref. [pQ|): 



ICI 

2 

K[ 

2 

+• 



lJu + l(() 2 -J»(()Ju + 2(()} + 



1 



ICIm 2 ^(Q 

2(z/ + l)(C 2 + / u 2 ) 

mom^o\ + t (c)J ?£V 



MO + ^+2(0 



MO 



2(f + l)|Cl 

^ +2 (C) 



J„(C)/ + ... 



(4.56) 



2(1/ + 1) V4(i/ + l)( I / + 2)|C| ICI 3 

Inserting this expansion into (|4.55[) , we find that the relevant integrals become analytically doable 
when v is taken large enough to make the integrals converge. In this way we have explicitly confirmed 
a number of terms in the above expansions of the massive spectral sum rules. 



5 Generalizations and outlook 



The small and large mass expansions for the partially quenched chiral condensate can of course in 
principle be extended to the two other major universality classes, corresponding to SU(2) gauge 
group with fermions in the fundamental representation, and SU(iV c > 2) with fermions in the adjoint 
representation. The challenge is here to find a convenient method that permits high-order expansions. 
In our present case all simplifications arose from the observation that the effective partition function 
is annihilated by two sets of Virasoro constraints, which in turn is rooted in the fact that this effective 
partition function is a r-function of an integrable KP hierarchy. The effective partition functions of 
the two other universality classes have not been nearly as well studied, although there are reasons to 
believe that they are what is known as "Pfaffian r-functions" (see, e.g., ref. |f2~lf ). It may thus be 
possible to derive analogous partition function constraints, which will permit high-order expansions 
with little labor. In this connection we also make the following observations. The Virasoro constraints 



( |2.S| ) and ( |3.41| ) both follow from the simple differential equation 

d^z^ . l 



BMikdMlj iJ4 

by changing variables to the pertinent expansion parameters. This differential equation by itself is 
a trivial consequence of the fact that the integration manifold is U(iVy + N v ), as it just amounts to 
inserting the unitarity condition U^U = 1 in the integrand of the effective partition function (2.1). 



Imposing two boundary conditions is thus sufficient to determine the partition function uniquely. In 
the case at hand it is the choice zjf* f \tf = 0) = 1 and Yq N * \tl = 0) = 1, together with the required 
property Zq N, '(M,M^) = Zq Nj \m^M) that provide the two conditions. 



Except for a change of the proportionality constant, here 1/4, the differential equation ( |5.ip also 
holds for the effective partition functions of the two other universality classes mentioned above, as 
well as for the one of QCD3. It is also instructive to notice that the restriction to a sector of fixed 
topological charge does not affect the differential equation. The boundary conditions imposed on the 
differential equation, however, separate the solutions, i.e., the partition functions. For instance, before 
the projection onto a fixed topological sector the partition function studied here depends on det(A4) 
and det(M^) as well as on MM.^ . 

Of course we can already at this stage obtain the partially quenched chiral condensates to lowest non- 
trivial orders for the two remaining ensembles by simply making use of the same small-mass expansions 
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of the partition function that were used to derive the first sum rules. This is a trivial exercise: from 
ref. H we learn that the effective partition function for SU(2) gauge group and Nf + N v fermions in 
the fundamental representation is 

ft = mM)Y (l + ^-^—^^M'M + ...), (5.2) 

where, now M. is an antisymmetric (2Nf + 2N V ) x (2Nf + 2N V ) matrix with the usual mass matrices 
+Ai and —M placed in the two off-diagonal blocks, and zeros in the two diagonal blocks. Using again 
eqs. ( 3.29 ) and (|3.30| ) we find the partially quenched chiral condensate to this order: 



^v{Pv) _ V Pv / r o\ 

E fi v 2(2N f + v - 1) + ' • ■ ' K ' 

which one can confirm matches the first term in the expansion of the result obtained from Random 
Matrix Theory (see the first of ref. [jllj). Note that this expansion has its first uncancelled de Wit-'t 
Hooft pole at 2Nf + v = 1, and the above term in the expansion is thus valid for 2Nf + v > 1. 

Similarly, for the universality class corresponding to Nf + N v adjoint fermions and arbitrary gauge 
group SXJ(N C ) we see from ref. Q that the partition function expansion is 

z > = AA{Mf f 1 + W7TJ^TT) TlM,M + • ■ ■) ' (M) 

where v = vN c is an integer. Using the present replica method this leads to a partially quenched 
chiral condensate of 

E fi v Nf + 2v + 1 

which one again can check matches the result from Random Matrix Theory [11]. There is no de Wit-'t 
Hooft pole in this case, and thus no retriction on the validity of this first term. 

Since one can obtain perturbative solutions to the partially quenched chiral condensate in all three 
universality classes, a natural question concerns the microscopic spectral density of the Dirac operator 
Ps{C> {w})- This density is given by the discontinuity of the partially quenched chiral condensate 
across the cut on the imaginary axis 0: 

Ps(C,{Vi}) = Disc| Mt _ if E„(/i„,{/ii}) = —Km[T, v (i( + e,{vi})-Z„(i(-e,{iJ,i})]. (5.6) 

(This identification holds when one considers T, u (fi v , {[J-i}) as a function of a real mass fj, v , and then 
replaces fj, v — * iQ ± e). It now seems straightforward to insert our small-mass and large-mass se- 
ries expansions for {p>i}), and derive corresponding series expansions for the spectral density 



ps(C; {pi})- But this is not possible. Let us first consider the small-mass expansion of eq. (|3.31|) . 
As we emphasized in section 3.1, this expansion is valid up to the first de Wit-'t Hooft pole. Except 
for the topological term v j \x v which always trivially yields the correct 5- function contribution to the 
microscopic spectral density, this power-series expansion of a finite number of terms has no cut across 
the imaginary axis. This may seem surprising, since the microscopic spectral density ps(C> {Pi}) does 
have a simple and well-defined power-series expansion. The explanation is, however, simple. Because 
of the de Wit-'t Hooft poles, we can derive the small-mass expansion of ^, u (fi v ,{fii}) up to (and in- 
cluding) order \x v f+2v . But the perturbatively expanded ps(C'>{Pi}) (see, e.g. the massless case 
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( pT27[) ) starts only at order £*N f +2u+l _ Thg 

microscopic spectral density is always precisely one step 
ahead of the order to which we can push the small-mass expansion of the partially quenched chiral 
condensate! In hindsight, a phenomenon like this had to occur, since the pure power series implied by 
the small-mass expansion cannot give rise to any discontinuity across the imaginary axis. So the de 
Wit-'t Hooft poles in fact precisely save what would otherwise be a paradoxical situation. 

The asymptotic large-mass expansion of section 3.2 cannot be used to extract the microscopic spectral 
density either. Here the reason is very different. Recall that the asymptotic expansion we have 
derived in section 3.2 is suitable \i v on the positive real axis. For example, for Nf = 1 our expansion 
simply coincides, up to an irrelevant normalization factor, with the asymptotic expansion of the Bessel 
function Io(m) : 

„„, <■-" ?(-')' r(* + i/2) 

This asymptotic expansion is correct, but it neglects exponentially suppressed terms of order / \fji 
and lower. Because of this, the asymptotic expansion ( pV7| ) does not reproduce the asymptotic expan- 
sion for Jo(/u) after rotation to the imaginary axis (the exponentially suppressed terms then become 
of magnitude comparable to those kept in ( |5.7| ), and that expansion is therefore no longer correct near 
the imaginary axis). We conclude that neither of the two expansions we have considered here are 
suitable for deriving the microscopic spectral density. 



6 Conclusions 

Using the Virasoro constraints on the effective partition function we have obtained small-mass and 
large-mass expansions for the QCD partition function valid in the scaling region V <C \jm\. In these 
expansions we can treat the number of valence quarks N v as a continuous parameter. This form is 
thus suited for a calculation of the partially quenched chiral condensate using the replica method. 
Two series of spectral sum rules for the QCD Dirac operator follow from the small-mass expansion. 
One series extends the Leutwyler-Smilga sum rules in a simple way to very high (here 8 th ) order, while 
the other series is new: it computes the spectral sum rules in the massive theory, and we have given 
series expansions in the physical masses for these new sum rules. In all cases we have checked, we 
have found complete agreement with earlier results based either on Random Matrix Theory or the 
super symmetric technique. The replica method constitutes a new and independent derivation of these 
results. 



With the expansions of the partition function at hand we have also derived the small-mass and large- 
mass expansions of the partially quenched chiral condensate. In all cases we have checked there is 
again complete agreement with earlier analytical predictions based on the two other methods. 

We have here restricted ourselves to small-mass and large-mass series solutions only because they 
provide the simplest framework in which to employ the replica method. This is not an inherent 
restriction, though. It is thus an open challenge to extend the method beyond these two series 
expansions. 
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